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A direct comparison of quantum and classical dynamical systems can be accomplished through the 
use of distribution functions. This is useful for both fundamental investigations such as the nature 
of the quantum-classical transition as well as for applications such as quantum feedback control. 
By affording a clear separation between kinematical and dynamical quantum effects, the Wigner 
distribution is particularly valuable in this regard. Here we discuss some consequences of the fact 
that when closed-system classical and quantum dynamics are treated in Gaussian approximation, 
they are in fact identical. Thus, it follows that several results in the so-called 'semiquantum' chaos 
actually arise from approximating the classical, and not the quantum dynamics. (Similarly, opposing 
claims of quantum suppression of chaos are also suspect.) As a simple byproduct of the analysis we 
are able to show how the Lyapunov exponent appears in the language of phase space distributions in 
a way that clearly underlines the difference between quantum and classical dynamical situations. We 
also informally discuss some aspects of approximations that go beyond the Gaussian approximation, 
such as the issue of when quantum nonlinear dynamical corrections become important compared to 
nonlinear classical corrections. 

PACS numbers: 03.65.-w, 05.45.Mt 



I. INTRODUCTION AND DISCUSSION 

A formal comparison of closed classical and quantum 
dynamical systems is conceptually difficult because what 
is natural in the classical description (trajectories) is not 
what is natural in the quantum description (amplitudes) 
and vice versa. While "semiclassical" limits of quantum 
mechanics are available (WKB, large N, coherent states, 
etc.) an important obstruction is the singular nature of 
the limit h ^ 0, and the fact that this limit docs not 
commute with the limit t —^ oo. Nowhere is this problem 
more important than in the study of quantum chaos. 

Here we adopt the viewpoint that the most natural 
language in which to compare closed classical and quan- 
tum systems is the language of phase space distribution 
functions as first set forth by Wigner jlj . Since quantum 
mechanics is a probabilistic theory such a description is 
physically attractive. Moreover, it has also turned out to 
be fruitful from a calculational perspective: Among sev- 
eral other applications, Wigner functions have been used 
to compute quantum corrections to classical statistical 
mechanics 0i Q' ^o clarify the WKB limit of quantum 
mechanics for both regular and stochastic orbits P, la, j 
to study quantum dynamics J}, and, in particular, dy- 
namical aspects of quantum chaos [3, Q including the 
effects of decoherence [l^, UM- However, we do wish to 
emphasize that results of this paper do not depend in 
any essential way on using the Wigner function: It is a 
convenient device to represent the density matrix, but 



* Electronic 



address: 



URL: http : //tSweb ■ lanl . gov/people/salman/ 1 



|habib@lanl.gov| 



certainly not necessary. We also emphasize that it is not 
our purpose here to discuss the quantum-classical tran- 
sition per se as would be the case in a measured (hence, 
open) quantum dynamical system |ia |. 

In the area of quantum chaos, the well-known work 
of Berry and Balazs Q is based on the idea of com- 
paring classical and quantum dynamics in phase space. 
These authors studied the evolution of a semiclassical 
(WKB) Wigner function corresponding to a classical 
curve in phase space, arguing that while at early times 
the Airy peak of this distribution function tracks the 
time-dependent classical curve, at later times, once the 
classical evolution produces structure on sub-7i scales in 
phase space, the quantum distribution is unable to con- 
tinue this tracking and there is a transition to a qual- 
itatively new regime. By identifying two sorts of non- 
linear structures, 'whorls' (generated by stable motion) 
and 'tendrils,' (generated by unstable motion) Berry and 
Balazs conjectured that this transition would occur on a 
timescale tc = 0{fi~^'^) for whorls and tc = 0(ln?i~^) 
for tendrils, the latter of which is the famous so-called 
'log' time. It is important to note that the existence 
of the log time does not necessarily imply the fail- 
ure of semiclassical methods beyond this time: later 
work has demonstrated that suitably improved semi- 
classical methods can be successful well past the log 
time [milllllll. 

A naive comparison of Wigner functions and classi- 
cal probability distributions is not possible as Wigner 
functions are not positive definite in general. This is a 
consequence of the fact that quantum mechanics admits 
interference effects, whereas classical mechanics does not. 
Nevertheless, one can set up a scheme to compare clas- 
sical and quantum dynamics in the following way: begin 



with an initial distribution function that is acceptable 
both classically and quantum mechanically (this condi- 
tion may be relaxed to allowing all acceptable quantum 
states, with or without a positive Wigner distribution) 
and then, (i) solve the two different classical and quan- 
tum Liouville equations, or equivalently, (ii) solve the 
corresponding infinite hierarchies of dynamical equations 
for the corresponding moments. The advantage of using 
the Wigner function is that, by introducing a c-number 
phase space description of quantum dynamics, one can 
directly compare the classical and quantum moments. 

A key feature of the Wigner formalism that we exploit 
here is the clear separation of kinematical and dynamical 
effects. Confusion has arisen in the literature as a conse- 
quence of the failure to recognize that H can occur in some 
places purely as a consequence of quantum initial condi- 
tions and, in other places, as a consequence of quantum 
dynamical corrections to classical dynamics. To sharpen 
one's arguments, especially in areas as controversial as 
quantum chaos, it is important to clearly distinguish be- 
tween these two situations. (Essentially this point, but 
phrased a little differently, has also been made by Ballen- 
tine, Yang, and Zibin [13 •) The importance of this last 
point is further underscored by the fact that classical and 
quantum mechanics have different underlying group sym- 
metries (symplectic vs. unitary) which accounts for the 
singular nature of the limit h ^f p!9 . Il9| . 

A clarification of just what is meant by such statements 
as "the classical limit" is important both conceptually 
and operationally (as emphasized in Ref. [i3|)- The usual 
textbook arguments involving wave packets and Ehren- 
fest's theorem attempt to compare a classical description 
in the sense of a single trajectory with a quantum ana- 
log, usually taken to be a sufficiently narrow wave packet. 
This is a rather limited point of view, especially in chaotic 
systems, where, due to their instability, the concept of in- 
dividual trajectories is not physically useful even classi- 
cally, and where quantum wave packet dynamics breaks 
down very quickly. It makes more sense, therefore, to 
compare classical distribution functions, describing an 
ensemble of trajectories with their appropriate quantum 
analogs as discussed above. (Such a viewpoint is in fact 
implicit in Feynman's discussion of the two-slit experi- 
ment 12(1 , a classic example of a highly nonlinear dynam- 
ical system.) The above philosophy of quantum-classical 
comparison has been implemented in some recent work 
on th e q uantum-classical corres pon dence in kicked sys- 
tems |2l|, the driven pendulum |23| . and the driven an- 
harmonic oscillator ^vf\. Such an approach JU| is also 
needed to analyze experiments in quantum chaos utiliz- 
ing ultracold atoms |23Ji in order to unambiguously dis- 
tinguish quantum from classical effects |2J| . 

Turning now to moment expansions |25| . the first im- 
portant point, discussed below, is that quantum dynam- 
ical effects do not arise in moments of less than third or- 
der. In an approach utilizing moment expansions around 
a centroid, this means that to first nontrivial order in fluc- 
tuations there are no dynamical quantum effects: This is 



the reason why the simple Gaussian approximation fails 
to sample any aspects of quantum dynamics. 

The infinite hierarchy of coupled equations for the mo- 
ments is almost impossible to deal with analytically; a 
variational approach is often attempted by assuming a 
particular functional form for the distribution functions, 
forcing thereby relations between lower and higher or- 
der moments (alternatively, one may utilize a cumulant 
expansion). A very popular approach is to assume the 
distribution function to be Gaussian which implies that 
all odd moments vanish and that all even moments may 
be expressed in terms of powers of the second order mo- 
ments. It is useful to distinguish between three differ- 
ent variants of the Gaussian approximation: The 'time 
dependent variational principle' (TDVP) method |26l |. 
the 'truncated Gaussian approximation' (TGA) of which 
Heller's method 27] is a non-self-consistent variant, and 
finally the multiple classical trajectory version of the sec- 
ond technique 13], which we will refer to as the 'multiple 
trajectory Gaussian approximation' (MTGA) (although 
in principle the MTGA can be applied to both types of 
single Gaussian approximations). To avoid confusion, the 
term 'Gaussian approximation' will be used when dis- 
cussing general properties of all three variants, otherwise 
they will be referred to specifically. 

It is easy to show that in TGA, and for the individ- 
ual packets in MTGA, no quantum dynamical corrections 
are included as only quadratic order moments are con- 
sidered. This is a generic problem with all Gaussian ap- 
proximations: in TDVP, while one does resum a certain 
series (which in a naive moment expansion is merely trun- 
cated) this resummation still does not incorporate terms 
arising from quantum dynamics. (As a consequence, it 
is not obvious whether the TDVP is any better than the 
truncated approximations.) For chaotic systems, both 
the TDVP and TGA will lose accuracy once the initial 
distribution becomes significantly non-Gaussian, which 
will happen at least on a Lyapunov timescale (typically 
this breakdown occurs much sooner) 28]. The MTGA 
is better behaved because a coherent sum of Gaussians 
can be a much better approximation for the full distri- 
bution function than a single Gaussian. The MTGA also 
keeps track of relative quantum phases between its con- 
stituent wave packets and this is the reason it can handle 
quantum revivals and other coherent effects. 

What is true for Gaussian approximations in nonlin- 
ear systems is also obviously true for linear systems where 
the Gaussian approximation is exact. One consequence of 
this is that the well-known independent oscillator model 
[29| , long considered a paradigm for quantum Brown- 
ian motion, in fact contains only classical dynamics. All 
quantum features in this model reside in the initial con- 
ditions. Once this key fact is realized, the derivation of 
the Master equation can be radically streamlined with 
no need to use path integral methods. Moreover, new 
light is shed on the derivation of classical Langevin equa- 
tions and the interpretation of precisely what is meant 
by 'noise' when the dynamics is not Markovian, yet the 



Master equation is local in time (A detailed discussion is 
given elsewhere |30|.') 

The main point of the present paper may be stated 
in the following way. Classically, or quantum mechan- 
ically, when one follows the evolution of a 'mean field' 
there are 'fluctuation' corrections (back-reaction) to the 
(classical) equations for the mean field. The quantum 
moment expansion described above allows one to sep- 
arate the back-reaction into its classical and quantum 
components. When the exact moment expansion is ap- 
proximated by either a quadratic truncation of the hier- 
archy, or by the Gaussian approximation, the quantum 
back-reaction is completely neglected. The only quantum 
feature still remaining is the uncertainty principle con- 
straint on the initial condition for the distribution func- 
tion. Therefore, if one wishes to compare the strengths 
of the quantum versus the classical back-reaction, the 
Gaussian approximation is wholly inadequate. These 
comments apply with equal force to Gaussian (Hartree) 
approximations which have attained a certain popular- 
ity in studying initial value problems in quantum field 
theory. Off-Gaussian corrections need to be taken into 
account, as is possible in principle, e.g., in a systematic 
1/A^ expansion oV\ (though this approximation has its 
own problems 32] ). We will discuss some of these as- 
pects further below. 

The Gaussian approximation has been applied to study 
chaos in quantum systems ^, .34] . In Ref . ,32] , the lead- 
ing order 1/A^ expansion (equivalent to a quadratic TGA 
in that particular case) was used to study a quantum me- 
chanical model which is the (0 -I- l)-dimensional limit of 
QED. Classical chaos was shown to exist in the dynam- 
ical equations for the correlation functions. In Ref. |34l |. 
the TDVP was applied to a one-dimensional quartic oscil- 
lator, and this system which is classically nonchaotic (in 
contrast to the situation in Ref. j33]), was shown, in the 
approximation, to display chaos. In interpreting these 
examples of "semiquantal" chaos, two questions imme- 
diately arise. The first is whether the chaos is just an 
artifact of the Gaussian approximation, and the second, 
whether it has anything to do with quantum dynamics. 
The existence of chaos in the general case of dynami- 
cal variation al ap proaches had been previously noted by 
Caurier et al [33 who argued that the onset of chaos sig- 
naled the breakdown of the approximation scheme. Mo- 
tivated by the results of Ref. 34] , the question was also 
tackled by Sundaram and Milonni j36i] using operator 
moments leading to a similar conclusion. In yet another 
investigation, it was also demonstrated that Gaussian ap- 
proximations break down before the timescale on which 
the approximate dynamics becomes chaotic J23. (For a 
recent study of semiquantum chaos in a different approx- 
imation, see Ref. 37].) The answer to the second ques- 
tion, given here, is that Gaussian semiquantum chaos, 
even if approximate, has nothing to do with quantum dy- 
namics. Thus the result of Pattanayak and Schieve |3J] 
is not counter-intuitive (how can "adding" quantum me- 
chanics to a nonchaotic classical system produce chaos?): 



it is a consequence of the approximation, and selective re- 
summation, of only the classical back-reaction terms in 
the moment expansion. This also rules out the heuristic 
explanation that their result is due to the inclusion some- 
how of quantum "noise," which then makes the system 
chaotic. 

The Gaussian approximation has also been applied by 
Zhang and Feng |3S] in an attempt to display an exam- 
ple of quantal suppression of chaos. Just as we reject the 
reality of semiquantal chaos, so must we reject such an 
interpretation of the calculation. What has been shown 
is that the classical back-reaction, in Gaussian approx- 
imation, can suppress chaos. This is not surprising in 
principle as the Gaussian approximation is completely 
uncontrolled and there is no guarantee that if the origi- 
nal system is nonintegrable, the Gaussian dynamics also 
will be (the reverse is true in semiquantum chaos where 
the original system can be integrable, and the Gaussian 
approximation, not). 

The above discussion leads naturally to another set 
of questions: What are the timescales on which non- 
Gaussian corrections become important? Are these 
timescales different for the classical and quantum back- 
reaction terms? As argued in Refs. \VA l3q it is cer- 
tainly reasonable to expect the TGA and TDVP to break 
down, for a classically chaotic system, on a timescale set 
by the largest Lyapunov exponent. Since ensemble av- 
eraged quantities do not possess chaotic signatures like 
a late-time Lyapunov exponent p39l ]. it is also clear that 
any time the Gaussian approximation is chaotic, the re- 
sult is simply wrong. The timescale set by this false 
chaos, in the cases where it occurs, is therefore precisely 
the timescale over which one expects the classical non- 
Gaussian corrections to become important. It is not a 
priori obvious what will happen for the quantum cor- 
rections, but given that the classical non-Gaussian cor- 
rections must eliminate the chaos on a (relatively) fast 
timescale, there is no dynamical imperative for these con- 
tributions to become dominant at early times. 

The concept of a crossover or 'break' time from clas- 
sical to quantum dynamics in closed systems (chaotic or 
not) is often invoked. However, the evidence for the util- 
ity of such a concept is surprisingly weak: The notion 
of the crossover time could well depend on the quanti- 
ties one chooses to compare, and on the dynamical ap- 
proximations one makes. If, for example, the estimates 
of this time are based on TGA or TDVP, then all one 
is computing is when the approximation breaks down. 
That this may not have anything to do with how well 
classical systems track quantum evolution (in a suitable 
coarse-grained sense), both for chaotic and regular sys- 
tems, has been pointed out in Ref. [I3]. Evidence for 
this point of view is contained, for example, in the work 
of Fox and Elston [2j on the kicked top and the kicked 
pendulum. These authors compare the time evolution 
of classical and quantum distributions, finding excellent 
agreement at early times (as expected, the Gaussian ap- 
proximation being accurate in this regime), but also at 



late times, when a coarse-grained quantum distribution 
is compared to the classical distribution. What this says 
is that the low order moments can be insensitive to the 
nature of the dynamics, despite the existence of nontriv- 
ial quantum dynamical effects. Investigations utilizing 
high-resolution simulations have shown that in some ex- 
amples of chaotic dynamical systems, while the classical 
and quantum moments do track each other very closely, 
there is a small 'mini-break' which occurs on a timescale 
longer than the dynamical time |l3- Whether this 'mini- 
break' has anything to do with log time arguments is not 
yet clear, though more recent numerical evidence appears 
to argue against this pO|. 



The set of equations of the simple, albeit non-self- 
consistent, TGA due to Heller [23 provide an instructive 
link to the computation of the actual classical Lyapunov 
exponents. Classical mechanics allows arbitrarily strong 
shrinking of a Gaussian ball in phase space around the 
fiducial trajectory; this shrinking is what allows the eff'ec- 
tive Hamiltonian for Gaussian fluctuations to reduce to 
that for linearized perturbations, which is just what de- 
fines the Lyapunov exponents. In this way, even though 
the classical Liouville equation is linear, the allowance of 
a singular trajectory limit permits the existence of chaos. 
However, in the quantum case, because of the finite value 
of h, the shrinking allowed classically is impossible be- 
yond a certain point, and the effective Hamiltonian for 
Gaussian fluctuations never becomes that for linearized 
perturbations. 



Finally, it is important to point out that the Gaussian 
approximation works extremely well in certain situations 
where localization of the quantum state is guaranteed, as 
can happen in studying the dynamics of observed quan- 
tum systems (quantum trajectories) (Cf. Ref. [ill)- In 
a related context, the use of a Gaussian ansatz for the 
Wigner function in quantum feedback control has also 
been shown to be very effective |4l|, |42, |43 , even when 
the true state can have substantial non-Gaussian fea- 
tures [43 . While this is not the focus of the present pa- 
per, some of the general formulae for Gaussian states dis- 
cussed here are also applicable to these, otherwise quite 
distinct, situations. 



The organization of this paper is as follows. In Sec- 
tion II, we study the classical and quantum Liouville 
equations, especially in the context of moment expan- 
sions. In Section III, we develop Gaussian dynamics us- 
ing moments and discuss various ramiflcations of the re- 
sults. We conclude in Section IV with a short summary. 



II. CLASSICAL AND QUANTUM DYNAMICS 
IN PHASE SPACE 



A. Liouville Equations 

Throughout this paper we will deal with continuous- 
time Hamiltonian dynamical systems, with Hamiltonian 



H 



P 
2m 



Vix), 



(1) 



where the potential V{x) is allowed to be explicitly time- 
dependent. We take the system to be one-dimensional: 
extension to higher dimensions is obvious. The classical 
equations of motion following from the Hamiltonian (Q 
are entirely equivalent to the classical Liouville equation 



dfc 
dt 



^cjc 



for the phase space distribution function fc{x,p), 
where the classical Liouville operator. 



m dx 



dVd_ 
dx dp' 



(2) 



and 



(3) 



Point trajectories can be recovered from the Liouville 
equation |(2Jl if one allows the distribution function to 
be a delta function over phase space. 

An important consequence of the classical evolution 
equation is the exact conservation of the quantities (when 
they exist; / is any phase-space distribution function) 



= / dxdp p 



(4) 



by the classical evolution equation (0), following from 
Liouville's theorem. (The trivial case n — 1 arises from 
the conservation of total probability.) This will acquire 
more signiflcance when we proceed to the quantum case. 
The corresponding quantum Liouville equation is writ- 
ten in terms of the Wigner function jj, which is a 'half- 
Fourier transform' of the density matrix p{xi,X2) in the 
coordinate representation. One introduces the sum and 
difference variables 



a; = (xi -fa;2)/2, A = xi - X2 
and then deflnes the Wigner function as 



fw{x,p) 



^/,(.,A)e-A/''dA. 



(5) 



(6) 



Details on the properties of this function may be found in 
Refs. 13,143. A particularly nice discussion, with special 
emphasis on conceptual pitfalls and the distinction be- 
tween kinematical and dynamical aspects, has been given 
by Tatarskii ^. (See also, Ref. gg.) 

For the purposes of this paper, it is sufficient to note 
the following: (i) fw is real and normalized to unity over 
the {x,p) phase space, (ii) fw is square integrable, with 



/ 



0-2 = / dxdp fw < TT^ 



2TTh' 



(7) 



the equality holding for pure states {27rfia2 — Trp^). 
This is a partial statement of the uncertainty principle 
in the Wigner formalism; localization of the distribution 
to a delta function on phase space, while classically al- 
lowed, is forbidden in quantum mechanics, (iii) The only 
pure state for which fw is positive everywhere is one with 
the Gaussian wave function (Hudson's theorem [4q) 



V' 



— aa; -\-bx-{-c 



(8) 



with a, 6, c complex and Re a > 0. (iv) Expectation val- 
ues of operators that are either pure powers in position or 
pure powers in momentum are directly given by the aver- 
ages of X™ or p™ taken with respect to fw- Expectation 
values of mixed operators like x^p^ cannot be obtained 
directly. The Wigner formalism is associated with Weyl's 
rule for the ordering of operators [43 : the phase space 
average of a mixed term like x"p'" corresponds to the 
expectation value of a completely symmetrized operator 
expression. 

The dynamics of the Wigner function is easily obtained 
from the quantum Liouville equation for the density ma- 
trix p{xi,X2)'- 

d 

ih,—-p{xi,X2) = [H(xi) - H{x2)*]p{xi,X2). (9) 

ot 
The Wigner transform of this equation gives 
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dx^ dp^ 



(10) 



where A 7^ 1, and with the assumption that the the po- 
tential can be Taylor-expanded. Often this is not the 
case and Eqn. IjlUI) becomes an integro-diffcrcntial equa- 
tion; we avoid this added complication in this paper. On 
comparing Eqns. H10|l and (01 we see that quantum and 
classical dynamics differ due to the last term in the quan- 
tum Liouville equation (|10|l . This apparent separation of 
quantum and classical pieces in the evolution equation is 
one of the original motivations for the Wigner function 
approach. However, since the definition of the Wigner 
function itself involves h, the actual situation is not as 
straightforward as it appears to be. In particular, the 
'quantum correction' in Eqn. H1U|I is not unitary. (Nei- 
ther is Lc, but their sum is. Ref. 18] contains a further 
discussion including the violation of positivity in the clas- 
sical approximation.) 

The quantum evolution equation H1U|) exactly con- 
serves the uncertainty principle constraint Q. This 
has a one-to-one correspondence with the classical evolu- 
tion preserving the corresponding classical quantity CT2 of 
Eqn. I0J. It is important to note that since the quantum 
evolution does not obey Liouville's theorem, in this case, 
&„ ^ 0,Vn ^ 1,2. Thus, since the conservation of cri 
is trivial in both classical and quantum evolutions, the 
constancy of only (72 in the quantum case acquires added 
significance. 



B. The Moment Hierarchies 

In the moment hierarchy approach to dynamics one 
attempts to replace a knowledge of the distribution func- 
tion by a knowledge of all-order values of the quanti- 
ties (x"), (p"), and (x"p™) where the bracket denotes 
an average taken over the distribution function. From 
now on ( )^ will denote an expectation value taken with 
respect to a classical distribution, and ( ) will denote 
an expectation value taken with respect to the quantum 
distribution function. Below we present an elementary 
derivation of the classical and quantum moment equa- 
tions using the evolution equations for the corresponding 
distributions (for alternative approaches, see Ref. 25] ). 

With the assumption that the distribution function 
and its derivatives are zero at the phase space boundary, 
in the classical case [using the classical Liouville equa- 
tion (0)] it is easy to show that : 



di 



(xV)^ = / dxdp a;V^ 
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(11) 



(12) 






(13) 



Of course, these equations are nothing but the Hamilton 

equations of motion averaged over the initial conditions. 

The quantum hierarchy is equally straightforward to 

obtain. The evolution of the spatial moments is given by 



formally the same as the classical case. However, 
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dt 
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where 



(-1)-^ 
e(n,A)=^ ' 



A! (n- A)! \2i 
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, (15) 



(16) 



differs from the corresponding classical equation. The 
upper limit of the sum in (|15|l is n, if n is odd, and 
n — 1 if n is even. Clearly the first dynamical quantum 
correction comes in only at third order, where 



d 
di 






Finally, the cross moments evolve according to 

k9fw 
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(18) 



As in Eqn. (|15|l . the upper limit of the sum is k if k is 
odd, and /c — 1, if it is even. 

A moment expansion around the centroid can now be 
developed by expanding in 'fluctuations' around 'mean 
flelds,' 



X 

P 









in the classical case, and 



{^), + s, 



p = {p)„ + v, 



(19) 
(20) 



(21) 
(22) 



in the quantum case. It is clear from the definition of the 
fluctuations that their expectation values vanish identi- 
cally. 

The strategy to follow now is simple: One begins with 
the lowest order classical dynamical equations from 111|) - 
(|13|l . or the corresponding quantum equations (|15l) - (|18|l . 
and substitutes in the appropriate pair from equations 
(|19|l - (|22|l given above. One then finds terms involving ex- 
pectation values of higher order fluctuations. Dynamical 
equations for these quantities can then be found by going 
to higher order in (^lj-(221) or (^Hl-llll)- Care has to be 
taken in interpreting the moment expansion: the classical 
and quantum mean fleld equations are not identical and 
the fluctuations in the quantum case contain both 'classi- 
cal' and 'quantum' contributions. An analogous moment 
expansion was performed by Sundaram and Milonni |3a | 
by working in a Heisenberg operator formalism. The mo- 
ment hierarchy described below is completely equivalent 
to the set of equations they derived. The advantage here 
is that we can compare directly with the corresponding 
classical hierarchy. 

To illustrate the above procedure we now carry out the 
moment expansion to quadratic order in the fluctuations. 
As will be seen below, to this order, there is no difference 
between the classical and quantum hierarchies (this is 
well-known in quantum optics, see, e.g., Ref. 50]) so it is 



sufRcient merely to present the classical situation. The 
equations for the average quantities are [set n = 1 in l|ll|) 
and ifT^ ]: 



-(x) = Mc 

dt ^ '^ m ' 

d IdV 



where 
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To obtain a closed set of equations we need to know the 
equations for the fluctuations correct to quadratic order. 
We begin by noting that 



which implies directly in turn 
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dt 



But we know already from the hierarchy that [set n — 2 
in HH)], 



dt 
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Comparing H27|) and 

d 
di 



in 



we find 
2 



{Sv)c 



(28) 



(29) 



which is the first of the required equations. Notice that 
the derivation of this equation is in fact exact as no terms 
were neglected. This is not true for the other two equa- 
tions to which we now pay attention. Proceeding as 
above, we first write 



then setting n = 2 in (|12|l . we find 



(30) 



l<^^>. 






= -2[{p).K^'^ + ^{p)An.v^'^ 
= ^(p)cjt{p)c-HvS),vP 



(31) 



Comparing H30() and (|31|l . we obtain 



dt 



{v')^^-2{vS)^V}'l 



(32) 



It is now apparent how to proceed to get the last equation 
and we merely write the final result 



|(.^).^^-(^^K^^^ 



(33) 



Equations ^, ^, ^, (EH, and 1(23) form a closed 
set of dynamical equations. Since at most only second or- 
der moments were used in the derivation, these five clas- 
sical equations are formally identical to their quantum 
counterparts. Had we truncated the moment expansion 
at one higher order there would have been a difference be- 
tween the quantum and classical equations since, as we 
have seen, the first quantum dynamical correction comes 
in only at third order. 

General expressions can be written for the exact clas- 
sical and quantum moment hierarchies by extending the 
procedure outlined above. In the classical case, restrict- 
ing attention only to those moments relevant to future 
discussion of the Gaussian approximation, one finds. 



!<*' 



>. = ^ («•'-'">, 



dt 



m 



c-)^(f). 



■q 



dx 



I <*"">' - 'M^ 



OX I m '^ 



(34) 



(35) 



(36) 



The quantum hierarchy is only slightly more compli- 
cated. For the same moments considered above. 



_(^«) = -(S''-^r]) 
dt ^ '1 m ^ '^1 



dt 
d 
dt 



(37) 



rii ^' '1 



'--'^^ 



.-idV\ 



+ A("\ 



dx / 



where 



n,n— 1 



A(") 



E Q("'^) 



A odd 






(38) 



(39) 



(40) 



Note that among the equations displayed, only H38|) 
contains quantum corrections compared to the classical 
equations. This is because H14II is already formally iden- 
tical to the corresponding classical equation, as is H18|l in 
the case fc = 1. 



III. THE GAUSSIAN APPROXIMATION 
A. Time Dependent Variational Principle 

We now proceed to the first variant of the Gaussian 
approximation. The idea here is to (i) assume that the 
distribution function is at all times a Gaussian (for more 
mathematical details on the formal properties of such 
states, see Refs. |5l|, Isll) and then, (ii) to write down 
the corresponding dynamical equations for the moment 
hierarchy. The approach given here is a generalization of 
the time dependent variational principle first considered 
by Dirac |2y| by considering mixed as well as pure states. 
The moment hierarchy approach is far more convenient 
than a direct application of a variational principle as car- 
ried out, for example, in Ref. 53]. As will be apparent in 
what follows, this approximation only samples the mo- 
ment hierarchy to quadratic order, and hence contains 
no dynamical quantum corrections. 

We assume the system to be always described by a gen- 
eral Gaussian distribution (here we make no distinction 
between classical and quantum) of the form 



/ = ^exp 



-a{x-x{t)f-b{p-p{t)f 
+c{x~x{t)){p~p{t))], (41) 



or in a more compact notation: 



/ 



1 



: exp 



j,k 



y/Det{2TTS) 

with xi = X, X2 ^ p, and the matrices 

„_i _ I 2a -c 
^ " ' -c 26 



K^i^J - ^3)SjkiXk - Xk) 



S - 



1 I 2b c 
Aab - c'^ \ c 2a 



(42) 
(43) 

(44) 



The distribution is normalized to unity over phase space, 



N = y/Det{2TTS) ^ 



27r 



VAab^ 



(45) 



and satisfies the constraint (jT)). This extra condition in 
turn leads to a constraint on the parameters of the dis- 
tribution: 



Aab 



(47ra2)' 



(46) 



This purely kinematical condition must be satisfied for 
quantum distributions. In general, classical distributions 
need not satisfy this condition, but if they do, then, as 
already discussed, this condition is exactly preserved by 
the classical dynamics. 

Given the distribution function specified above, it is 
trivial to show that 



(x) = X, (p) = p, 
and that the quadratic moments are 



{S') - ^11 



1 
1 



(77(5) =: 512 = 



c 
46 



^-'--a 



An^a^ 



2^2 



47rV| 



(47) 

(48) 
(49) 
(50) 



4a/3 (47rcr2)2' 
The kinematical constraint H46|) can now be expressed as 

1 



{rj'){6') = {vSr + 



(4^^2)2' 



(51) 



which is a statement of the uncertainty principle gener- 
alized to mixed states. For pure states, (72 = l/iirh, and 
(|51|l takes the familiar form (the wave function is now 
that of some Gaussian state, more general than just a 
squeezed coherent state) 



{v'){s') = {vsf + '-^. 



(52) 



All odd-order moments are zero and the higher order 
even moments can be computed from the quadratic ones 
by a simple application of the result, valid for Gaussian 
distributions, 



{{Xj 



,M^)).Y.''^J?M\, 



\ dxk 



(53) 



where x = (xi,a:2) and ^(x) is some arbitrary function 
of X. The moments relevant to us are 



(<52») = 



(2n-l)! 



j.m\ 



2"-i(n 

(2n-l)! ^ 2.„ 
2n-i(n-l)! V' / ' 

(2n+l)! 



2\" 



2"n! 



{V6){S') 



(54) 
(55) 
(56) 



where ()51|l has been used to simplify the last expression. 
The above expressions are the only ones needed to obtain 
the moment hierarchy in Gaussian approximation. 

The dynamical equations under the Gaussian approx- 
imation are easy to obtain. First, 



dt 



M 

m 



(57) 



and 



d d_ 

dt ^P^ = dfP 



dt 

dV 
dx 



P_ 
m 



(58) 



2! 



y(3) 



3! 



y(4) 



y^ \li_y(2n+l) 



n!2'i 



(59) 



where we have used (|54() in the last step. The dynamical 
equation for (^^^ jg nothing but (|29|l which is already 
exact. The equation for (jy^^ follows from ()35(l or H38|l : 



I*"') - 






= -2 



r7iy(i)+<5V^(2)^^^(3) 



3! 



°° T/(2fe+2) t. 



fc=0 



2'^-ifc! 



(60) 



This result has been written in a simple form by using 
(|56|l along with the constraint H51() . The last equation 
we have to deal with is for {5rj), this following from 136|l : 



dt 



(Sv) 



1 

m 



W) 



dx 



- (ry2) - A f T/(i) + 5^(2) + ^^3) 
m ^ ' \ \ 2! 



3! 



n+l 



1 /„2\ _ y- lii ^(2«+2) 



(61) 



Note that we have used (|^ in the last step. 

We now have five equations, two for the mean fields, 
and three for the fluctuations. It is easy to check that 
the constraint l(5^ is exactly preserved by the dynami- 
cal equations. Therefore the constraint may be used to 
eliminate one of the three dynamical equations for the 
fluctuations. For example, if we choose to keep (^S'^'^ and 
{Sj]) as the preferred variables, then the two dynami- 
cal equations for the fluctuations are (|29(l and 1)61(1 with 
(ry2^ written in terms of (6^) and (Srj) using 1)5 l|l . Notice 
that these equations were obtained without having to go 
to moments higher than those of order two. Therefore 
they do not contain any quantum dynamical contribu- 
tions, just as advertized. Of course, h can still appear in 
these equations, but in a purely kinematical sense. This 



will happen if the initial state is pure (cr2 — l/27r?i) and 
we use the constraint to eliminate one of the dynamical 
equations for the fluctuations. 

To conclude our discussion of the TDVP, we note that 
this approximation is conservative and further, can be 
mapped to a classical Hamiltonian for a new set of dy- 
namical variables [^. The first point can be easily 
checked by writing the expectation value of the Hamilto- 
nian (here assumed explicitly time independent) in Gaus- 
sian approximation, 



(^> = ^ (/) + (n^)) 



2m 



Y(2k) 



ipf + W)\+Y.-^{n (62) 



/c=0 



2'=^fc! 



and then verifying, using the dynamical equations given 
above, that 



>- 



(63) 



It is interesting to note that any consistent truncation 
of the TDVP is also conservative (by consistent trunca- 
tion, we mean: keep only the moments up to some finite 
order, i.e., go up to the same finite order in all the equa- 
tions). As it happens, Heller's approximation 27| is not 
consistent in this sense because it does not account for 
the back-reaction on the mean field thereby missing a 
quadratic term that should have been kept; consequently 
it fails to be conservative. 

It is easy to show that the four dynamical equations for 
the centroid and two chosen fluctuation variables are in 
fact Hamiltonian. Eliminating (rf^ from H61|) using the 
constraint H51|l and introducing the new variables J54] 34] 



{5') 



2\ — 



{ri5) = p-i, 



(64) 



in H62|) . we find that the Hamiltonian 



Hr, = (H) = 



1 2 1 2 

°° „2n 

P TA(2n) 



1 



2mp'^{ATTa2) 



E^^y^'^'i^) 



(65) 



generates the correct equations of motion as can be easily 
verified. It follows trivially from (|65|l that any consistent 
TGA is also Hamiltonian. 

The expectation value of the Hamiltonian (|62|l is pre- 
cisely the quantity that is minimized in the Gaussian 
effective potential approach 55] . Since ?i appears only in 
kinematical guise in 1)62(1 it follows that no essential quan- 
tum effects (i.e., those without formal classical analogs) 
can be seen in this approximation. For example, tunnel- 
ing/overlap corrections in a double-well system are not 
accessible in simple Gaussian approximation {56| . 



B. Semiquantum Chaos and Related Issues 

In Ref. [sj, Pattanayak and Schieve studied a one- 
dimensional anharmonic oscillator in Gaussian approxi- 
mation. (See also Ref. /citcC2.) Even though the origi- 
nal system is clearly not chaotic, in the approximation it 
can be, because we now have four coupled nonlinear dy- 
namical equations, rather than the original two. This is 
the main technical point of their paper and is, of course, 
correct. The authors stressed the surprising nature of 
their result, that dynamics of expectation values may be 
chaotic, "even though the system has regular classical 
behavior and the quantum behavior has been assumed 
regular" jsj- Their conclusion has been criticized by 
Sundaram and Milonni [Sg who point out that the Gaus- 
sian approximation becomes unreliable very quickly when 
either the classical system is chaotic, or when the Gaus- 
sian dynamical equations are themselves chaotic. The 
more general theoretical setting is the following: Any 
variational approximation to the dynamics of a quan- 
tum system based on the Dirac action principle leads to 
a classical Hamiltonian dynamics for the variational pa- 
rameters [53 ■ Since this Hamiltonian is generically non- 
linear and nonintegrable, the dynamics thus generated 
can be chaotic, in distinction to the exact quantum evo- 
lution |28l l35| ; the Gaussian approximation represents a 
particularly simple special case. 

Once we realize that, in fact, the analysis presented 
in Refs. [33 and |3J| is purely classical, the results are 
easy to understand. First, they are not surprising in the 
sense that quantum mechanics is not somehow catalyz- 
ing chaos. Second, we see why the approximation breaks 
down as it does: if the original system is chaotic then, 
'tendrilization,' in the sense of Berry and Balazs [g will 
happen on a timescale set by the maximal Lyapunov ex- 
ponent by which time the Gaussian approximation is def- 
initely wrong. If, on the other hand, the original classical 
system is not chaotic, but the Gaussian set of equations 
are, all this tells us is that the non-Gaussian corrections 
are becoming very important on the timescale set by the 
maximal Lyapunov exponent of the false chaos, and need 
to be included to obtain the correct classical dynamics. 
Third, the explanation that quantum noise is somehow 
promoting the chaos is not tenable since the Gaussian 
approximation contains only classical fluctuations. 

An ostensibly opposite result was obtained by Zhang 
and Feng j38j : These authors implemented the TDVP ap- 
proach for the kicked rotor, finding a suppression of the 
classical chaos. This, they claimed, was a simple exam- 
ple of a generic suppression of classical chaos by quantum 
mechanics. Our interpretation of their result is different: 
since it is a consequence of merely classical dynamics, 
what they have really displayed is a suppression of chaos 
by the classical back-reaction, computed in Gaussian ap- 
proximation. Thus, there really is no information about 
quantum suppression (or even of a classical 'fluctuation 
suppression,' since the calculation is only approximate). 
Real quantum suppression of classical chaos in this model 
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is discussed in Refs. 1581 15 



C. The Truncated Gaussian Approximation 

The full Gaussian approximation (TDVP) involves 
summations to all orders of the fluctuations (though 
these are somewhat trivial as the higher terms that are 
summed are disconnected pieces, i.e., no cumulants of 
order higher than two are encountered). This approx- 
imation differs from a truncation approach where one 
merely truncates the moment hierarchy at some finite 
order. We can incorporate the truncation idea in Gaus- 
sian approximation by requiring the distribution to be 
Gaussian but still stopping the summations at some finite 
order. As discussed above, this truncation is conserva- 
tive and Hamiltonian in the sense that the corresponding 
equations of motion for dynamical variables can be de- 
rived from a Hamiltonian. The TGA may be useful in 
some situations because of its computational simplicity. 

An example of this idea is the approach due to 
Heller [23. In our language, Heller's instructions are 
the following: (i) Drop the fluctuation corrections to 
the equation for the average momentum (|59|) . and (ii) 
keep only the leading order (quadratic) corrections in the 
equations for the moments [Eqns. l|^ . (|^ . and l|CT|) ]. 
The motivation for this approach as originally stated is 
to provide a semiclassical limit that does not quite cor- 
respond to the usual h ^ limit of WKB theory (For 
example, the quadratic order TGA is exact for the har- 
monic oscillator, but WKB is not). As will become clear 
in what follows, dynamically this is still a classical ap- 
proach and its ostensible quantum features are easily ex- 
plained: (i) the appearance of h is purely kinematical, 
and (ii) "penetration" of the Gaussian into classically for- 
bidden territory is simply a classical consequence of the 
distribution function having support at high momenta. It 
is no surprise that this approximation has trouble with 
tunneling (which is accessible in a WKB approach) since 
it does not incorporate quantum dynamical corrections. 

By carrying through with the prescription given above, 
we obtain the dynamical equations: 



d 

di 
d_ 

dt 
d_ 
dt 



dt 
d _ 
Jt^ 



to' 



.yd), 

- {H 






(66) 
(67) 
(68) 
(69) 
(70) 



The first two equations 162|l and l|63() are just the clas- 
sical equations of motion for a particle in the potential 
V{x). This approximation thus ignores the back- reaction 
of fluctuations on the mean field. As for the fluctuations 



themselves, they are sensitive only to the second deriva- 
tive of the potential, as in the quadratic approximation 
of Section IIB. Heller's method is not equivalent to a 
quadratic truncation because there is no back-reaction 
on the mean flelds and because it still enforces the con- 
straint (|51|l . reducing the number of dynamical equations 
from flve to four. However, it is not a consistent quadratic 
TGA because of the failure to include the back-reaction. 
Just as in the TDVP and quadratic approximations, this 
approximation only includes classical dynamics. 

As a consequence of the lack of self-consistency at 
quadratic order, the expectation value of the energy is 
not conserved: 



4(^) = My(3)(^2^ 



dt 



2m 



(71) 



It is a trivial matter to improve this approximation by 
keeping the back-reaction consistently to quadratic order. 
(Note, however, that this approximation is conservative 
for linear systems.) Even though this method is not uni- 
tary, it has been applied very successfully, especially in its 
multiple trajectory variant, to chaotic systems p^ll4l|. 

The implementation of the TGA in the Schrodinger 
picture also involves taking into account a time vary- 
ing global (Maslov) phase for the wave function. Since 
the Wigner function is a transform of the density ma- 
trix, it does not contain this phase information. If there 
was more than one packet (as in MTGA), however, the 
Wigner function would keep relative phase information 
for which there is no classical analog. 



D. Connection with Lyapunov Exponents 

A straightforward observation is that the set of equa- 
tions H66|) - H70|l are exactly those needed to determine 
the classical Lyapunov exponent, with {x,p) specifying 
a flducial classical trajectory, and the fluctuation vari- 
ables specifying the linearized evolution of a Gaussian 
ball of initial conditions around the fiducial trajectory. 
As before, imposition of the constraint arising from the 
constancy of (72 leaves only two equations for the fluctu- 
ations, and these are Hamiltonian separately, with 



Hf = ^1^ + 
2m 



1 



2mp^(4...)^ + ^^'^(^)^^- (^^) 



This would be precisely the linearized Hamiltonian for 
the system were it not for the second term which en- 
forces the constraint. The definition of the Lyapunov 
exponent requires the double limit of first squeezing the 
initial Gaussian ball onto the fiducial trajectory (in all 
phase space directions) and then taking the limit t ^> 00. 
The first limit corresponds to taking CT2 — > 00 (clearly 
not allowed quantum mechanically for ti finite), which is 
the same as dropping the second term in Hp. Inciden- 
tally, this is also a good way to see how chaos can occur 
in the Liouville theory despite the fact that the classical 
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transport equation is linear (like the Schrodinger equa- 
tion or the quantum Liouville equation): it arises simply 
because the delta function or trajectory limit is allowed 
classically, but not quantum mechanically. 



E. The Multiple Trajectory Gaussian 
Approximation 

The idea behind the MTGA is to improve the TGA 
by including multiple classical trajectories. One does 
this by propagating more than one Gaussian, using the 
TGA for each individual packet. There are two reasons 
why this is a good idea. First, a sum of Gaussians can 
track the classical distribution even at times when the 
distribution is highly non-Gaussian (although this may 
involve a high computational cost as very many Gaus- 
sians might be needed in the sum). The second issue 
involves quantum coherence in phase space. The classi- 
cal version of this approximation simply sums over the 
individual Gaussians to give the composite distribution 
function. This is not true in the quantum case: a co- 
herent sum of Gaussian wave functions, it turns out, is 
not a sum of the associated Wigner functions, but rather 
a sum augmented by "off-diagonal" quantum coherence 
terms 'jj\ . (These off-diagonal terms account for the Airy 
fringes in the WKB approximation for the Wigner func- 
tion discussed by Berry (J.) Thus the MTGA contains 
nontrivial quantum information, and is in this sense, un- 
like TDVP and TGA, not a classical approximation. 



IV. SUMMARY 

The comparative evolution of c-number classical and 
quantum moments yields an obvious method to compare 



classical and quantum mechanics as dynamical theories. 
In doing so here, we have restricted attention largely 
to the case of the Gaussian approximation. The time- 
dependent Gaussian approximation in quantum mechan- 
ics has a direct, completely equivalent, classical analog 
in terms of an ensemble of classical particles also treated 
in the same approximation. Thus every 'quantum' re- 
sult in this approximation has an exact classical coun- 
terpart. Part of the reason for the misinterpretation of 
'semiquantum' chaos as a quantum dynamical effect was 
simply due to a too naive comparison with classical me- 
chanics: a comparison against trajectories, rather than 
against swarms of trajectories. Using the more general 
approach presented here, it becomes clear that, in this 
context, h plays no fundamental role. The key point 
is that a dimensionful constant also appears in classical 
transport theory: a2, as defined in Eqn. Q, which plays 
the role of an effective h there - such a constant is miss- 
ing from the equations of motion of classical trajectories 
but arises once classical distributions have to be evolved 
by the Liouville equation. 
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